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Abstract 

Instrumentation and methods for rapid screening and selection of plants with improved water use efficiency are 
essential to address current issues of global food and fuel security. A new imaging system that combines chlorophyll 
fluorescence and thermal imaging has been developed to generate images of assimilation rate (A), stomatal con- 
ductance (g s ), and intrinsic water use efficiency (WUE,) from whole plants or leaves under controlled environmental 
conditions. This is the first demonstration of the production of images of WUE; and the first to determine images of 
g s from themography at the whole-plant scale. Data are presented illustrating the use of this system for rapidly and 
non-destructively screening plants for alterations in WUE, by comparing Arabidopsis thaliana mutants (OST1-1) that 
have altered WUE, driven by open stomata, with wild-type plants. This novel instrument not only provides the poten- 
tial to monitor multiple plants simultaneously, but enables intra- and interspecies variation to be taken into account 
both spatially and temporally. The ability to measure A, g s , and WUE; progressively was developed to facilitate and 
encourage the development of new dynamic protocols. Images illustrating the instrument's dynamic capabilities are 
demonstrated by analysing plant responses to changing photosynthetic photon flux density (PPFD). Applications of 
this system will augment the research community's need for novel screening methods to identify rapidly novel lines, 
cultivars, or species with improved A and WUE; in order to meet the current demands on modern agriculture and 
food production. 

Keywords: Chlorophyll fluorescence imaging, dynamic responses, leaf heterogeneity, screening, thermal imaging, water use 
efficiency. 



Introduction 

One of the greatest challenges plant scientists currently face 
is global food and fuel security. Water availability is a major 
constraint of crop yield (Sinclair and Rufty, 2012) and is the 
single most important factor limiting food production, with 
significant yield losses reported under water deficit (Boyer, 
1982; Mueller et ol, 2012; van Ittersum et al, 2012). In order 



to combat the predicted impacts of increasing drought epi- 
sodes on crop yield, there is an urgency to identify plants 
and the underlying mechanisms for improved water use effi- 
ciency (WUE). Unfortunately, a major constraint for such 
crop improvements is the ability to monitor WUE rapidly 
and non-destructively. It is essential that new techniques 



Abbreviations: A, C0 2 assimilation; C a , extracellular C0 2 concentration; Q, intracellular C0 2 concentration; F', steady-state fluorescence measured under light; 
F m ', maximum fluorescence measured under light; F q 7F m ', photosystem II operating efficiency; g c , cuticle conductance; g b leaf conductance to water vapour; 
g B , stomatal conductance; GUI, graphical user interface; / g , index of stomatal conductance; IRGA, infrared gas analysis; /WUEj, imaged intrinsic water use efficiency; 
PPFD, photosynthetic flux density; r aw , leaf boundary layer resistance to water vapour; r HR , parallel resistance to heat and radiative transfer; r m , leaf resistance to 
water vapour; s, slope of the curve relating saturation vapour pressure to temperature; T Dry , mean temperature of the dry standard; 7~ LeElf , mean temperature of the 
leaf surface; r Wet , mean temperature of the wet standard; VPD, vapour pressure deficit, WUE, water use efficiency; WUEj, intrinsic WUE; y, psychrometric constant; 
6, temperature phrase normalizing leaf temperature with the wet and dry standards. 
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and approaches are developed for phenotyping (Fiorani and 
Schurr, 2013) and to screen for limitations in WUE. 

Traditional methods of determining WUE that quantify 
plant yield or biomass relative to the amount of water used 
(Stanhill, 1986; Bacon, 2004; Passioura, 2004) are unsuitable 
for rapid screening for several reasons. These agronomic tech- 
niques are not only destructive but also rely on an integrated 
measurement of biomass/yield at the end of the growing 
season relative to the amount of water used over the growth 
period (Chaerle et al, 2005; Morison et al, 2008). Carbon 
isotope discrimination (Farquhar et al, 1982) has been suc- 
cessfully used to identify crop cultivars with greater WUE 
(Condon et al , 2004); however, this technique also relies on 
an integrated measure of WUE over a period of plant growth. 
Additionally, the technique does not provide an indication 
of whether differences in WUE are driven by C0 2 assimila- 
tion (A) or water loss (Farquhar et al, 1989; Jones, 2004ft), 
although the incorporation of oxygen isotope measurements 
can provide an indication of rates of evaporation from the leaf 
surface (Farquhar et al, 1998; Barbour, 2007). Leaf-level gas 
exchange measurements of the rate of A relative to transpira- 
tion provide an immediate and non-destructive measure of 
instantaneous WUE (Penman and Schofield, 1951) or 'intrin- 
sic water use efficiency' (WUE ; ), when stomatal conductance 
(g s ) is used instead of transpiration as a measure of water 
loss (Meidner and Mansfield, 1968; Jones, 1998). Although 
this approach is flexible in term of the time scale of when 
measurements can be made, an infrared gas analyser (IRGA) 
can only take singular measurements on one plant or leaf, at 
one point in time. Thus, all of the techniques described above 
to assess WUE have limitations as screening tools, as they 
tend to be time-consuming and/or destructive. Additionally, 
biomass and carbon isotope measurements only provide a 
lifetime measure of WUE based on cumulative seasonal con- 
ditions, which may mask specific phenotypic traits that could 
be advantageous in future breeding programmes (Weyers and 
Meidner, 1990; Chaerle et al, 2005). Using a combined chlo- 
rophyll fluorescence and thermography imaging approach, a 
non-invasive, high-throughput, high resolution tool has been 
developed to screen WUE ; based on calculated images of A 
and g s produced from measurements of photosynthetic effi- 
ciency (F q '/F m ') and leaf temperature, respectively. 

Chlorophyll fluorescence has long been used to examine 
various photosynthetic parameters in leaves (Baker, 2008) 
and it is well established that the operating efficiency of pho- 
tosystem II (PSII; F q '/F m ') is related to changes in C0 2 assim- 
ilation in leaves (Genty et al, 1989, 1990; Krall and Edwards, 
1990; Cornic and Ghashghaie 1991; Edwards and Baker, 
1993; Siebke et al, 1997). However, the relationship is com- 
plex and depends on the surrounding gaseous environmental 
conditions. In C 3 plants, the relationship between the operat- 
ing efficiency of PSII is linearly related to photosynthetic C0 2 
fixation of leaves, but only when photorespiration is inhibited 
and C0 2 assimilation represents the major sink for the end- 
products of electron transport, namely ATP and NADPH 
(Baker and Oxborough, 2004; Baker, 2008). Chlorophyll fluo- 
rescence images taken on C 3 plants under a low [0 2 ] (20 mmol 
moL 1 ) can be converted to images of C0 2 assimilation using 



suitable calibrations (Di Marco et al, 1990; Genty et al, 
1990; Cornic and Ghashghaie, 1991; Cornic, 1994; Genty and 
Meyer, 1995). 

Infrared thermography (IRT) provides a powerful imag- 
ing tool for rapidly, non-invasively, and remotely measur- 
ing leaf temperature as a surrogate for g s (Omasa et al, 
1981; Hashimoto et al, 1984; Jones, 1999). Leaf tempera- 
ture depends on evaporative cooling, and is a function of g s 
(Jones, 1992). For example, leaf temperature increases as sto- 
mata close and restrict evaporative water loss. Thermography 
has been used for rapid screening of genotypic variation in g s 
(Wang et al , 2004) as well as early diagnosis of drought stress- 
induced changes in g s (Grant et al, 2006; Morison et al, 
2008). In the last 20 years there has been increasing interest in 
quantitative evaluation of stomatal conductance from meas- 
urements of leaf temperature, using the basic energy balance 
equations (Jones, 1999, 2004a; Leinonen et al, 2006; Grant 
et al, 2007). The use of thermography to determine g s has 
been optimized through the development of standard proto- 
cols which take into account the surrounding environment, 
and even the distribution of stomata between the two leaf 
surfaces (Jones, 1999; Guilioni et al, 2008). Thermography 
has become a standard technique to determine g s in both 
glasshouse (Grant et al , 2006) and field environments (Grant 
et al, 2007). 

Here the development of a novel imaging system that incor- 
porates measurements of chlorophyll fluorescence and ther- 
mal imaging under controlled gaseous conditions is described. 
Chlorophyll fluorescence images of F q '/F m ' provide a quanti- 
tative image of A, whilst images of leaf temperature are con- 
verted to g s using well-defined methods that take into account 
the surrounding environment. Further manipulation of these 
two images provides for the first time an image of intrinsic 
water use efficiency (IWVE^A/g^). Previous researchers have 
used combined chlorophyll fluorescence and thermal imaging 
approaches to evaluate photosynthetic performance in relation 
to stomatal behaviour (e.g. Chaerle et al , 2005, 2007; Lawson, 
2009; Glenn, 2012), but the majority of these studies have been 
carried out at the leaf or tissue scale (Omasa and Takayama, 
2003; Messinger et al, 2006) and have not been used to deter- 
mine WUE, but have been focused on physiological analysis 
of mechanisms that coordinate responses observed between 
mesophyll photosynthesis and stomatal behaviour. A major 
advantage of the combined imaging approach reported here 
is the ability for multiple samples to be measured at any one 
time and the fact that spatial heterogeneity within plants and 
leaves can be readily identified. Differences between wild- 
type (WT) and open stomatal mutant (OST1-1; Merlot et al, 
2002) Arabidopsis thaliana plants are demonstrated using 
the system, with spatial and temporal heterogeneity in A, g s , 
and TWUEj being observed in the images. It is important to 
take into account such heterogeneity as it is well established 
that g s and photosynthesis are not uniform over a leaf sur- 
face (Oxborough and Baker, 1997; Weyers and Lawson, 1997; 
Weyers et al, 1997; Mott and Buckley, 1998; Lawson et al, 
2002; Peak et al, 2004; West et al, 2005; Kamakura et al, 
2012) and that such heterogeneity is also dynamic (Lawson 
and Weyers, 1999) often being driven by variations in the 
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microenvironment (Flexas and Medrano, 2002; Peak et al, 
2004; Lawson et al, 2012). The system has not only been 
constructed to image spatial differences in WUE ; , but it has 
been specifically designed to facilitate dynamic measurements, 
which allow the impact of changing environmental conditions 
on stomatal behaviour to be assessed in relation to photosyn- 
thetic performance and WUE. 

Materials and methods 

Plant material 

Arabidopsis thaliana genotypes Columbia-0 (Col-0), Wassilewskija-0 
(Ws-0), and Landsberg erecta (Ler), and the mutant Open Stomata I 
(OST1-1) were grown in a controlled environment at 23 °C and 1.1 
kPa vapour pressure deficit (VPD) day and night. The photoperiod 
was 8/16h light/dark with a photosynthetically active photon flux 
density (PPFD) of 1 35 ± 10 umol nT 1 s _1 . Two-week-old seedlings were 
transferred either to 100 cm 3 pots or to 6-well culture plates (Nunc, 
Roskilde, Denmark) containing compost (Levington's F2S, Everris, 
Ipswich, UK). Plants were maintained under well-watered conditions. 
A layer of vermiculite (Vermiculite Lite, Sincair, UK) was placed over 
the compost surface of plants grown in well plates to improve the con- 
trast between the plants and the background during imaging. 

Phaseolus vulgaris L. cv. 'Evergreen' were grown in a temperature- 
controlled glasshouse at 23 ± 4 °C. Lighting was supplemented by 
sodium vapour lamps (600 W; Hortilux Schreder, The Netherlands) 
when external solar radiation fell below 500 umol m 2 s PPFD, 
during a 10 h period. Plants were grown from seed in 650cm 3 pots 
containing compost (Levington's F2S) and watered every 2 d with 
Hoagland's nutrient solution. 

Combined chlorophyll fluorescence and thermal imaging system 

A chlorophyll a fluorescence imaging system (Fluorlmager, 
Technologica, Colchester, Essex, UK), previously described by 
Barbagallo et al. (2003), was modified by repositioning the camera 
from being directly above the chamber to a 90 ° angle, whilst main- 
taining the same distance from the subject material (Fig. 1). A silver- 
coated mirror (Thor-Optics, Dachau, Germany) was hinged on an 



axis directly above the original camera port. At a 45 ° angle, the mir- 
ror reflected the chlorophyll a fluorescence signal directly onto the 
camera. For thermal imaging, a thermal camera (TH7100 Thermal 
Tracer, NEC Avio Infra-red Technologies Co. Ltd, Japan) was posi- 
tioned in the original location of the chlorophyll fluorescence cam- 
era, directly above the imaging port. Pivoting the mirror allowed 
thermal and chlorophyll fluorescence images to be captured within 
2 s of each other (Fig. 1). The thermal camera has a temperature 
resolution of 0.1 °C, and all measurements were made at a distance 
of 0.45 m and emissivity (e) of 0.98. 



Gas control in imaging chamber 

In order to control concentrations of 0 2 , C0 2 , and H 2 0 vapour 
during imaging, an in-house designed chamber was constructed. As 
the spectral wavelength used for chlorophyll fluorescence excitation 
and thermal and fluorescence emission span from the visible to the 
infrared (450-14 000 nm), a chamber window could not be used and 
instead an open-top design was employed. The chamber was built 
from Perspex with inner dimensions of 145mm (length)x 105mm 
(width)x95mm (depth), including a 10mm width flange on the top 
surface. With the exception of the base, the chamber consisted of 
an inner and outer wall separated by a 10 mm gap. The outer walls 
were connected on each of the four sides by 6 mm PTFE tubing con- 
nections. The inner wall was perforated with 1 mm diameter holes 
at a density of 9 per 100 mm 2 , which was optimal for maintaining 
homogenous gas concentrations whilst minimizing leaf movement 
through turbulence. Total gas flow entering the chamber was typi- 
cally 0.87 1 s -1 . Within the chamber, target gas concentrations of N 2 , 
0 2 , and C0 2 were individually maintained by mass flow control- 
lers (EL Flow, Bronkhorst, Ruurlo, The Netherlands), connected 
to compressed gas cylinders containing 100% N 2 , 0 2 , and C0 2 , 
respectively (British Oxygen Company-Industrial Gases, Ipswich, 
UK). In order to control water vapour concentration, gas was bub- 
bled through temperature-regulated gas wash bottles (Cole-Parmer, 
London, UK) prior to entering the chamber. 

Gas composition in the chamber was monitored at plant height by 
sampling air with a diaphragm pump (Type 124, ADC Hoddesdon, 
Herts, UK) at 500 cm 3 min . Oxygen concentration was measured 
with a flow-through oxygen sensor (S101, Qubit Systems, Kingston, 
Canada) which was calibrated using 0 2 -free air and a 205 mmol 
mob 1 [OJ standard (British Oxygen Company-Industrial Gases). 




Fig. 1 . Schematic diagram of the system used to image whole plants for chlorophyll a fluorescence and temperature under controlled 
conditions. The imaging system was modified to allow the attachment of two cameras, the thermal camera being directly positioned 
above the plant and the fluorescence camera situated at 90 ° to the thermal camera, utilizing a silver-backed mirror at 45 ° to capture 
images of F q 7F m '. The plant was positioned within an open-topped chamber (photograph insert) where concentrations of 0 2 , C0 2 , and 
H 2 0 were maintained at a typical flow rate of -0.87 I s~ 1 N 2 , with all gases passing through a humidifying system prior to entering the 
chamber. Concentrations of 0 2 , C0 2 , and H 2 0 were measured every second by an IRGA and oxygen sensor. 
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Both C0 2 and H 2 0 vapour concentrations were measured with an 
IRGA (Li- 840, Li-Cor, NE, USA) calibrated weekly using a stand- 
ard gas for C0 2 (±2.5% tolerance) (British Oxygen Company) and a 
dewpoint generator (LI-610; Li-Cor) for H 2 0 vapour. 

Since F q /F m ' is only linearly related to A when photorespira- 
tion is inhibited, it was essential that [OJ within the chamber could 
be reduced and maintained at 20 mmol moL 1 during the imaging 
process and immediately returned to ambient concentration when 
complete. Figure 2a shows that [OJ could be rapidly reduced from 
205 mmol mol -1 to 20 mmol moL 1 within <50 s of switching the 
gas input. In this example, low [OJ were maintained (within 2 mmol 
moL 1 of the target value) for several minutes, although typically 
<20 s were required to capture an image of F'IF m ', Once an image 
was taken, the [0 2 ] was returned to ambient within 30 s (Fig. 2a). 



Figure 2b illustrates the stability of [COJ (400 umol m 



±5%) 



and H 2 0 concentrations (±10% of target values) during the [OJ 
changes and for the duration of the experimental procedure. 



Estimating carbon assimilation from chlorophyll fluorescence 
parameters 

F q '/F m ' is calculated from measurements of steady-state fluorescence 
in the light (F') and maximum fluorescence in the light (F m ') since 
F q '/F m '=(F m '-F')/F m ' . Images of F' were taken when fluorescence 
was stable at the desired PPFD, whilst images of maximum fluores- 
cence were obtained after a saturating 800 ms pulse of 5500 umol 
m -2 s -i ppFD (Oxborough and Baker, 1997, 2000; Baker et al. , 2001). 
The intensity of the saturating pulse is sufficient to saturate PSII in 



relatively low PPFD-grown plants, but may need to be increased or 
modified for high PPFD-grown samples (see Loriaux et al, 2013). 
Using the first fully expanded leaves of 5-week-old A. thaliana (Ws-0), 
F'/F m ' and A were measured at 15 different C0 2 concentrations. 
Leaf intercellular [COJ (Cj) and A were measured using an IRGA 
(CIRAS-1, PP Systems, Amesbury, USA) and F q IF m ' was deter- 
mined from simultaneous images of chlorophyll fluorescence. These 
data were used to produce A/Q and F q '/F m '/Q response curves, that 
were used to calibrate F q '/F m ' with net C0 2 assimilation (Morison 
et al., 2005). The standard cuvette window of the IRGA was 
replaced with non-reflective glass to enable fluorescence images to 
be taken of leaves. At each external C0 2 concentration (C a ), Q and 
A were allowed to reach steady state before images were captured. 
Measurements started at ambient C a of 400 umol mol -1 , before C a 
was decreased step-wise to a lowest concentration of 50 umol mol -1 , 
then increased step-wise to an upper concentration of 2000 umol 
mol -1 . Leaf temperature and VPD were maintained at 25 °C and 1.2 
kPa, respectively. Four replicate A/Q and F q /F m '/Q response curves 
were measured at 200, 500, and 800 umol m -2 s -1 PPFD. Plots of 
F q 'IF m ' against A over the Q range and at these different PPFDs 
were used to determine the relationships between measured A and 
F q IF m . 



Estimating stomatal conductance from leaf temperature 

IRT has been widely recognized to be a powerful imaging tool, 
capable of rapidly and non-invasively measuring leaf temperatures 
(Jones, 1999, 2004a; Leinonen et al, 2006; Grant et al, 2007). As 
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Fig. 2. Concentrations of (a) 0 2 and (b) C0 2 (filled circles) and water vapour (triangles) in the measuring chamber during an experiment. 
Oxygen concentration was switched from atmospheric (210 mmol mol" 1 ) to 20 mmol mol -1 five times while maintaining C0 2 and H 2 0 
vapour concentrations. 
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leaf temperature is dependent on evaporative cooling, it can be used 
as an indirect measure of leaf conductance to water vapour (gj) or 
its reciprocal, leaf resistance (r w ) (Jones, 1992, 1999; Guilioni et al, 
2008). Although this measure of leaf conductance also includes 
the water lost through the cuticle, termed cuticle conductance (g c ), 
the value is relatively small and usually neglected, subsequently the 
terms leaf and stomatal conductance (g s ) are used interchangeably 
(see Jones, 1999). Leaf temperature also depends on the environmen- 
tal conditions around the leaf (Jones, 1992, 2004a) and, in order to 
determine g s from thermal images, an estimation of the boundary 
layer resistance to water vapour (r aw ) is needed along with known wet 
and dry temperature reference standards (Jones, 1999). Two tempera- 
ture references are required, one that provides an infinite resistance 
to water vapour (e.g. leaf material greased on both sides), whilst the 
second provides a near-zero resistance to water vapour (e.g. leaf sur- 
face painted with a detergent-water mix). The temperature standards 
normalize the measured leaf temperature to the environmental con- 
ditions surrounding it and it is assumed that these surfaces have the 
same radiative properties (Jones, 2004a). The application of the wet 
standard to either one or both sides of the leaf is used to account for 
the distribution of the stomata (Guilioni et al, 2008). In this study, a 
one-sided wet standard was used for Phaseolus vulgaris leaves and a 
two-sided wet standard was used for A. thaliana leaves. 

Leaf resistance to water vapour (r w ) was calculated from thermal 
images (Guilioni et al. , 2008) for anisolateral leaves using the follow- 
ing equation: 
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where r aw is the leaf boundary layer resistance to water vapour, 
s is the slope of the curve relating saturation vapour pressure to 
temperature, and y is the psychrometric constant. r HR represents the par- 
allel resistance to heat and radiative transfer, and 8 is the temperature 
phrase, normalizing leaf temperature with the wet and dry standards; 
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where 7i e af> 7w et , and T Dly are the mean temperatures of the leaf wet 
and dry standards, respectively. This temperature phrase normalizes 
leaf temperature (T Le!l {) between a surface of lower resistance and 
temperature ( r Wet ) and a surface of higher temperature and a resist- 
ance greater than that of the leaf (7D ry ). Estimates of g s were made 
from the reciprocals of r w ; g s =ll r w . 

In order to test how robust the method was for estimating g s from 
leaf temperature, IRGA measurements of g s were taken alongside 
estimates determined from images of leaf temperature from P. vul- 
garis and A. thaliana (WS-0 and Col-0). f aw values were calculated 
using damp filter paper leaf replicates with areas of 0.0025 m 2 and 
0.00044 m 2 for P. vulgaris and A. thaliana, respectively, following the 
vapour-flux density method of Weyers and Meidner (1990). 

The rate of water loss from the leaf replicates was determined 
under the imaging system from the change in weight (±1 mg), meas- 
ured every 30 s for a period of 15min. The average surface tem- 
perature of the filter paper was analysed using thermal imager 
software (Radiometric Thermography Studio Complete, Metrum, 
Wokingham, UK) to monitor temperature stability throughout the 
measurement period. VPD was maintained at 1.7 kPa (±0.006 SE) 
and air temperature at 21.7 °C (±0.04 °C SE). Concurrent measure- 
ments of air speed (average 0.11 m s~') over the filter paper were 
made using a hot-wire anemometer (Model 425, Testo, Alton, 
Hampshire, UK). 

Leaves of P. vulgaris were kept flat by a wire support, with the 
adaxial surface facing upwards. Five-week-old A. thaliana plants 
were maintained in pots large enough to ensure leaves were parallel 
to the soil surface. A 30% Tween (Sigma, Gillingham, Kent, UK) 



solution or vacuum grease (Dow Corning, Midland, Michigan, 
USA) were applied to different 1 cm 2 areas of the leaf to create the 
wet and dry standards, respectively (see above). To generate a range 
of stomatal conductances, plants were left for 30min to reach steady 
state at PPFDs ranging from 100 umol m~ 2 s -1 to 1000 umol nr 2 s . 
The plant was then placed under the imaging system and a thermal 
image was taken followed immediately by an IRGA measurement. 
Conditions inside the imager were maintained at 1.7 kPa VPD and 
21.7 °C air temperature. 

Construction of images of water use efficiency 

The rationale and approach to constructing images of WUE ; 
(/WUEj) is illustrated in Fig. 3. Image values of F q '/F m ' (Fig. 3a) 
and leaf temperature (Fig. 3b) are converted to A (Fig. 3c) and g s 
(Fig. 3d) using the calibrations described above. To produce the 
image of JWUE;, the A image (Fig. 3c) was rotated, scaled, and 
interpolated to re-map spatially pixel values of A to values of g s . At 
the individual pixel level, values of A were divided by values of g s to 
produce pixel values of /WUE; (Fig. 2e; A/g s =rWUE{). 

Software development for image construction 

In-house specifically designed software (ImFluTem; IFT) with 
a graphical user interface (GUI) was developed using Matlab 
(Mathworks, Natick, MA, USA) to facilitate data processing and 
construction of images of WUE;. Leaf temperature and F q 'IF m ' 
data were downloaded from their respective commercial image 
capture programs and stored in a x,y,z matrix format which were 
uploaded into the Matlab program as text files. The maximum 
image area for leaf temperature images was 150 cm 2 (240 x 320 pix- 
els) whilst chlorophyll fluorescence images were smaller (90 cm 2 ) 
with a greater pixel resolution (700 x 520 pixels). The program 
allowed users to select the PPFD used and to map the data to 
A using the appropriate calibration (see Fig. 4). It also allowed 
manual input of other environmental variables (including air tem- 
perature, wet and dry reference temperature, and boundary layer 
resistance) required for the conversion of leaf temperature data into 
g s . Alternatively, an option was available that allowed user selection 
of areas of the images that corresponded to wet and dry standard 
reference material. As mentioned above, A images were rotated, 
scaled, and interpolated within the program to align pixels between 
the two images. Although the majority of the scaling and rotation 
values were embedded in the coding of the program, manual input 
options were designed to allow users full control over the mapping 
of images. Frequency distribution of the raw data provided a ref- 
erence check for the range of interpolated and calculated values 
of the new images produced, with the removal of any pixel values 
falling out of the original distribution. Images of /WUE ; were con- 
structed from individual pixels values of A divided by their cor- 
responding pixel value of calculated g s to produce pixel values of 
/WUE;. Values were mapped to a colour pallet that corresponded 
to 10 evenly distributed data bins determined by user-defined maxi- 
mum and minimum values. Output data were stored as text files 
along with corresponding JPEG images of A, g s , and /WUE;. An 
additional GUI assisted with processing multiple images (particu- 
larly imported when dynamic protocols were employed) and area 
selection for analysis. Multiple files were uploaded using a wildcard 
input function. The images and data generated used the same num- 
bered wild card function as the original file. 



Screening for differences in water use efficiency 

Three (4-week-old) plants of A. thaliana WT (Ler) and three 
OST1-1 mutants grown in a well plate were dark adapted for 20min 
before being placed in the imaging system and exposed to a PPFD 
of 200 umol nr 2 s -1 . After a minimum of 20min and at stable F' , the 
[0 2 ] was decreased from 205 mmol moL 1 to 20 mmol moL 1 for ~50 s 
and a thermal image recorded. This was immediately followed by the 
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Fig. 3. Typical (a) F q 7F m ' and (b) temperature images of an A. thaliana plant illustrating the production of (c) C0 2 assimilation (A), 
(d) stomatal conductance (g s ), and (e) intrinsic WUE (A/VUE|) images. Calibrations were applied to convert pixel values within each image 
into /A and g s data, respectively. The wet and dry temperature standards used in the calculation of g s are indicated by arrows on the 
temperature image. To produce an image of WUEj (/A/g s =AA/UE i ), the A image was rotated, scaled, and interpolated onto the image 
of g s . The colour bar beneath each image shows the range of parameter values. For further details of this process, see Materials and 
methods. 



application of a saturating pulse to produce F m " allowing calcula- 
tion of F q '/F m '. Following these measurements, [0 2 ] was returned 
to 205 mmol mol -1 . A second set of measurements was taken after a 
further 15min. All images were made at a [C0 2 ] of 400 umol mol -1 , 
VPD of 1 . 1 kPa, and air temperature of 2 1 °C. A 1 cm 2 greased circle 
of leaf material provided the dry standard (r Dry ) while damp filter 
paper, the same shape and size as a fully expanded leaf, provided the 
wet standard (7w et ). Estimates of /WUE; were made as described 
above. 



of Q values (35-1400 umol mol -1 ). Air temperature and VPD were 
maintained at 25 °C and 1 .2 kPa, respectively. Immediately after an 
IRGA reading was taken, the plant was rapidly transferred to the 
imaging system that was maintained at the same environmental con- 
ditions as the leaves in the IRGA chamber, with the exception of 
[0 2 ] that was maintained at 20.5 mmol mol -1 . When F' was stable, 
thermal and fluorescence images were taken and used to calculate 

/WUE;. 



Monitoring of dynamic changes in water use efficiency 

To image the dynamic response of A, g s , and /WUE ; to a step- wise 
increase in light, 5-week-old A. thaliana plants were positioned in 
the imaging chamber at 200 umol m 2 s PPFD, 400 umol mol -1 
[COJ, 1.1 kPa VPD, and 21 °C air temperature. A 1cm 2 patch of 
grease was applied to the abaxial and adaxial surface of a leaf to 
provide the dry standard, whilst damp filter paper was used for the 
wet standard. Every 3 min, [0 2 ] was reduced from 205 mmol mor 1 
to 20.5 mmol mol -1 and thermal and fluorescence images taken. 
After 15 min, PPFD was increased to 800 umol nr 2 s and images 
recorded every 3 min for a further 30min. To compare A, g s , and 
/WUEj values calculated from the images with those obtained using 
standard IRGA, measurements were made concurrently on similar 
plants of the same age grown and treated under the same conditions. 

Comparison of water use efficiency determined by imaging and 
gas exchange 

IRGA measures of WUE; were made on 5-week-old A. thaliana 
plants at 200 umol nr 2 s" 1 or 800 umol nr 2 s -1 PPFD and a range 



Results 

Relationship between Fq'/Fm' and C0 2 assimilation for 
leaves in imaging chamber 

The response of F q '/F m ' and A to changing internal [C0 2 ] (Q) 
at 20 mmol moF 1 [0 2 ] and at the different PPFDs is shown 
in Fig. 4a and b. As photorespiration was suppressed, both 
A and F q '/F m ' showed similar shaped saturation functions of 
Cj, with an initial linear increase as Q increases before they 
plateau at a given Q, which is dependent on PPFD. A plot 
of A against F q '/F m ' (Fig. 4c) shows a robust linear relation- 
ship between the two parameters at each of the three PPFDs. 
Regression analysis between the two parameters provided the 
functional relationship with which measurements of F q 'IF m ' 
were converted to A, and the regression coefficients indi- 
cated that >94% of the variation was accounted for in the 
relationships. 
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Fig. 4. The response of C0 2 assimilation (a) and PSII operating efficiency (b), estimated from F q 7F m ', to changes in internal C0 2 
concentration (Q) under 20 mmol mor 1 0 2 and at 200 (open circles), 500 (grey circles), and 800 (filled circles) ^mol rrr 2 s~ 1 PPFD. The 
correlation between C0 2 assimilation and PSII operating efficiency (c) was fitted using a linear regression for each PPFD intensity. All 
measurements were made on 5-week-old A. thaliana. Air temperature and VPD were 25 °C and 1 .2 kPa, respectively. Data are the 
means with standard errors (n=3-5). 
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Validity of estimating stomatal conductance from 
thermal images 

Stomatal conductance calculated from images of leaf tem- 
perature were compared with independent measurements of 
g s obtained by infrared gas exchange analysis (Fig. 5). These 
data demonstrated a significant correlation between meas- 
ured and calculated g s ranging from 40-640 mmol m 2 s -1 to 
640 mmol m 2 s -1 , for both plant species. This relationship 
validates the use of thermography to evaluate g s accurately. 

Demonstration of use of fluorescence and thermal 
imaging to detect differences in water use efficiency 

Images of A (Fig. 6a, d) and g s (Fig. 6b, e) determined from 
images of F q '/F m ' and leaf temperature, respectively, were 
used to calculate images of TWUE ; (Fig. 6c, f) for WT and 
OST1-1 A. thaliana measured after 20min and 35min of 
applying 200 umol nr 2 s" 1 PPFD. After 20min (T 2 q) in the 
light, there was no significant difference in A between the WT 
and OST1-1 (Fig. 6a), with all plants showing a mean value 
of 3.1 umol m 2 s" 1 . OST1-1 exhibited a lower A than the WT, 
although the difference was not significant. The distribution 
of pixel values within each image indicated that the majority 
of OST1-1 values were below 4 umol m 2 s -1 (Supplementary 
Fig. SI available at JXB online), whereas WT values covered 
a greater distribution range (0.5-8 umol nr 2 s -1 ) and were 
more evenly distributed within it. Stomatal conductance was 
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Measured Stomatal Conductance (mmol m 2 s" 1 ) 

Fig. 5. A comparison between stomatal conductance calculated 
from thermal images with that measured using an IRGA. To 
stimulate a range of conductances, leaves of Phaseolus vulgaris 
(filled circles) and A. thaliana thaliana Col-0 (grey circles) and WS-0 
(triangles) were exposed to PPFDs between 200 pmol rrr 2 s~ 1 and 
2000 jimol rrr 2 s _1 . The solid line represents a 1 :1 relationship 
(P < 0.0001 , fl s =0.96). 



significantly (T>=0.069) higher in OST1-1 plants (Fig. 6b); 
mean g s values for OST1-1 were 395 mmol nr 2 s _I which were 
31% greater than the average WT values. This difference was 
also reflected in the distribution of values within the images 
(see Supplementary Fig. SI). The mutant showed a frequency 
distribution that was skewed toward higher conductance val- 
ues (> 500 mmol nr 2 s -1 ) compared with the WT, where the 
greatest frequency of pixel values were <200 mmol nr 2 s -1 . 
The higher g s value observed in the OST1-1 plants, in con- 
junction with no difference in A, resulted in a significantly 
(T>=0.019) lower TWITE; compared with the WT. The average 
7WUE, in WT plants was 44% greater than in the mutants 
(Fig. 6c). After a further 1 5 min (T 35 ) at 200 umol nr 2 s" 1 PPFD, 
A in both plant types increased by 1.9 to 5.2 umol nr 2 s" 1 
and 4.8 umol nr 2 s" 1 for WT and OST1-1 plants, respectively. 
Average g s in the WT had increased by 43% ( 1 1 8 mmol nr 2 s -1 ); 
however, there was no change in mean OST1-1 g s (Fig. 6e). 
The higher g s and A after 3 5 min in WT plants (Fig. 6d) 
resulted in these plants having a similar TWUEj to the OST1-1 
mutants (Fig. 6f). The frequency distribution of pixels 
reflected the considerable variation observed in the images, 
with values ranging from 0.5 umol rrr 2 s" 1 to 10 umol rrr 2 s _I 
for A, 500 mmol nr 2 s _1 to 1000 mmol nr 2 s _I for g s , and 
0.0025 umol C0 2 /mmol H 2 0 nr 2 s" 1 to 0.05 umol C0 2 /mmol 
H 2 0 nr 2 s -1 for TWUEi (Supplementary Fig. SI). However 
no difference in distributions were observed between the two 
plant types at 35 min. 

Imaging dynamic changes in water use efficiency 

WUE; can oscillate rapidly with changes in both A and g s , 
driven by fluctuations in the environment and in particular 
light. To demonstrate the effect of changing light on stomatal 
behaviour and the impacts on A and WUE ; , an A. thaliana 
plant was subjected to a step-wise increase of 600 umol nr 2 s" 1 
PPFD following stabilization at 200 umol nr 2 s" 1 PPFD. In 
order to verify that the values obtained from the images were 
typical of those obtained using standard gas exchange meth- 
ods, simultaneous measurements of A, g s , and WUE; were 
captured using an IRGA (Fig. 7) and the combined imaging 
system (Fig. 8). There were no significant differences between 
measurements obtained using the two methods. A increased 
from 6.8 umol nr 2 s" 1 to 11.9 umol nr 2 s -1 at 800 umol nr 2 s _I 
PPFD (Fig. 7a). At 200 umol nr 2 s 1 PPFD, steady-state 
mean g s was 165 mmol nr 2 s _1 and 196 mmol nr 2 s _I in imaged 
and IRGA measured plants, respectively. Stomatal conduct- 
ance increased ~30% in both measured and calculated data 
after 30 min at 800 umol nr 2 s" 1 PPFD (Fig. 7b). Although 
the response of g s determined from the imaging system mir- 
rored that of IRGA measurements, the values on average 
were 13% lower. This discrepancy is most probably due to 
differences in the size of the boundary layer surrounding 
the leaves under the imaging system and inside the IRGA 
cuvette (see Discussion). Steady-state WUE ; increased from 
-0.036 umol C0 2 mmol/H 2 0 nr 2 s" 1 at 200 umol nr 2 s" 1 PPFD 
to 0.052 umol C0 2 /mmol H 2 0 nr 2 s" 1 at 800 umol nr 2 s" 1 
PPFD in measurements obtained using the IRGA (Fig. 7c). 
TWUEj was slightly higher than IRGA measurements at 
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Fig. 6. Screening of mutant OST1 -1 (left column of three plants) and wild-type Ler (right column of three plants) for differences in C0 2 
assimilation (/\), stomatal conductance (g s ), and A/VUE,. The plants were dark adapted for 20min before the PPFD was increased to 
200 jimol rrr 2 s~ 1 . Images were taken after 20min (a-c) and after 35min (d-f). Images of A (a and d), gs (b and e), and /WUE, (c and f) 
were calculated. All images were captured at 400 jimol mor 1 [C0 2 ], 21 °C air temperature, and 1 .1 kPa VPD. The colour bar between 
each image shows the range of parameter values. 



200 ]xmo\ m~ 2 s _1 PPFD, with an average value of 0.043 umol 
C0 2 /mmol H 2 0 m~ 2 s _1 . /WUE; steadily increased after 
PPFD was increased to 800 umol nr 2 s" 1 PPFD and reached 
a maximum value of 0.072 umol C0 2 /mmol H 2 0 m -2 s -1 
before stabilizing at a slightly lower value of 0.056 umol C0 2 / 
mmol H 2 0 m~ 2 s" 1 , identical to the IRGA measurements. 
/WUEi was significantly greater (7'=0.02) at 800 umol m~ 2 
s" 1 due primarily to the significant and rapid increase in A 
and the slower, smaller increase in g s following the step-wise 
increase in PPFD (Fig. 7c). The small non-significant differ- 
ences observed between measured and calculated values were 
most probably due to the heterogeneity observed between 
leaves and across plants. This variation cannot be taken into 
account by IRGA measurements of individual leaves. Images 
taken at 9, 21, and 30min illustrate considerable variation in 
A, g s , and /WUEj within and between leaves of an individual 
plant at any one time point (Fig. 8), as well as the change in 
parameter values with time and PPFD. At 200 umol m~ 2 s -1 
PPFD, ~70% of A values were between 5 umol nr 2 s -1 and 
9 umol nr 2 s _1 . When PPFD was increased to 800 umol nr 2 
s -1 , the median A value increased from 6.1 umol nT 2 s -1 to 
10.3 umol nT 2 s _I along with increased variation in the dis- 
tribution of pixel values (1-30 umol nT 2 s -1 ; Supplementary 
Fig. S2a at JXB online). The variation in A and g s resulted in 



heterogeneous patterns of /WUEj; for example, older leaves 
generally exhibited a 25% lower 7WUE; under 800 umol nr 2 
s _1 than younger leaves (Fig. 9; Supplementary Fig. S2). 

Comparison of water use efficiency determined by 
imaging and gas exchange 

A direct comparison between WUE; determined from the 
combined imaging system and measurements taken using 
an IRGA (Fig. 9) showed a strong positive correlation. The 
measurements were taken under different environmental con- 
ditions including a range of PPFDs (200-800 umol nT 2 s -1 ) 
and a range of Q values (35-1400 umol mok 1 ) on different 
leaves from different A. thalicina plants. Measurements of 
WUEj varied substantially with these different conditions, 
with values ranging from 0.008 umol C0 2 /mmol H 2 0 nT 2 s~' 
to 0.07 umol C0 2 /mmol H 2 0 nT 2 s -1 . 

Discussion 

A new imaging system capable of near-instantaneous com- 
bined chlorophyll fluorescence and thermal imaging has been 
developed in order to generate images of A, g s , and TWUEj 
from attached individual leaves and whole plants under 
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Fig. 7. Changes in measured (filled circles) and calculated (open circles) values of C0 2 assimilation, stomatal conductance, and 
WUE/AA/UEj were determined on leaves of 5-week-old A. thaliana plants during a stepwise increase in PPFD from 200 jimol rrr 2 s~ 1 to 
800 ^mol rrr 2 s~ 1 PPFD after 15 min (f). The asterisk (*) denotes data points presented as images in Fig. 9. Air temperature and VPD 
were 25 °C and 1 .2 kPa, respectively. Data are means with standard errors (n=3). 



strictly controlled environmental conditions. This is the first 
demonstration of the production of images of WUEj and pro- 
vides the prospect of rapid and non-destructive screening of 
plants for alterations in A, g s , and WUEj. Using this system, 
the impact of stomatal behaviour on 7WUE! was examined 
by comparing WT and known stomatal mutant A. thaliana 
plants and by monitoring plant responses to changing PPFD, 
as a driver of stomatal behaviour. 

Values of g s indicating the extent of stomatal open- 
ing (Omasa et al, 1981; Jones, 1992) have previously been 
inferred from thermal images (Omasa et al, 1981; Croxdale 
and Omasa, 1990; Inoue, 1990; Taconet et al, 1995; Jones, 
1999; Costa et al, 2013); however, to date, only one study 
has produced images of g s based on thermography. Omasa 
and Takayama (2003) analysed abscisic acid (ABA)-driven 



spatiotemporal changes in g s in conjunction with fluorescence 
measurements of F q '/F m ' and non-photochemical quenching 
(NPQ) using a combined imaging system. The authors pro- 
duced images of g s from measurements of leaf temperature; 
however, these images only examined a small area of a sin- 
gle leaf (30 x 30mm). This present study also illustrated the 
potential for quantitative analysis of spatial and temporal 
variation in g s and fluorescence in intact leaves that could pro- 
vide valuable information regarding the coordination of sto- 
matal and photosynthetic responses (Omasa and Takayama, 
2003; West et al, 2005; Aldea et al, 2006; Chaerle et al, 
2007). Values of g s estimated from thermal images have been 
found to be closely correlated with direct measurements of g s 
from both porometry (Jones, 1999) and IRGA (Fig. 5) meas- 
urements, consequently giving a high degree of confidence in 
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Fig. 8. Images of C0 2 assimilation, stomatal conductance, and /WUE, taken at 9, 21 , and 30min as shown in Fig. 7. Arrows indicate 
older leaves selected for analysis. Air temperature and VPD were 25 °C and 1 .2 kPa, respectively. The colour bars on the right show the 
range of parameter values. 



the images of g s generated from this system. This has enabled 
spatial and temporal variation in g s to be routinely moni- 
tored, allowing heterogeneity in /WUE; to be assessed when 
g s images are mapped to images of A (Fig. 3). 

In 'unstressed' A. thaliana plants, spatial variation in A, 
g s , and TWUEj was apparent within and between measure- 
ments of individual leaves, even under stable gas concentra- 
tions, PPFD, and temperature (Fig. 6, 8). However, once 
an intentional perturbation is introduced, such as the step- 
wise increase in PPFD, this variation in /WUE; was greatly 
increased. These differences in magnitude and rate of change 
were mostly driven by g s . It is well established that stomatal 
responses are an order of magnitude slower than A, which 
often results in a disconnection between A and g s following an 
alteration in the environment (Lawson et al, 201 1, 2012), that 
manifests itself in spatial and temporal variation in A and g s 
(Barradas and Jones, 1996; Weyers et al, 1997; Lawson and 
Weyers, 1999; Lawson et al, 2002). The degree of variation 
and the patterns observed are not surprising and are simi- 
lar to many previous studies that have reported spatial and 
temporal variation in either A or g s or other related variables 
(reviewed by Pospisilova and Santrficek, 1994; Weyers and 
Lawson, 1997; Lawson and Weyers, 1999) at scales that range 
from leaf to canopy (Weyers et al, 1997) and in response to 
various abiotic and biotic stresses (Tang et al, 2006; Ehlert 
and Hincha, 2008; Scholes and Rolfe, 2009; Bauriegel et al, 
2011; Nabity et al, 2012). A lack of coordination between A 
and g s and a lag in stomatal behaviour of between 5 min and 
lOmin is observed when g s does not initially change as PPFD 
is increased, whereas A responds immediately to increasing 



PPFD (Fig. 7). The response of A generally occurs within 
lmin, while the response of g s only occurs after a lag of 
several minutes and may take tens of minutes to complete 
(Meidner and Mansfield, 1968; Barradas and Jones, 1996; 
Wilmer and Fricker, 1996; Lawson and Weyers, 1999). It 
should however be noted that this is not the only explana- 
tion for variation in photosynthetic capacity (Miranda et al , 
1981); leaf anatomy (Sharkey, 1985; Terashima, 1992), leaf 
temperature (Hashimoto et al, 1984), boundary layer thick- 
ness (van Gardingen and Grace, 1991), and water relations 
(Slavik, 1963) all possibly play a role. Such spatial and tem- 
poral heterogeneity accentuates the value and benefit of using 
an imaging approach to A, g s and WUEj over alternative tra- 
ditional cuvette-based methods, which are generally confined 
to taking a single reading on an individual leaf or area of 
leaf which may or may not represent the average of the entire 
plant. The close relationship between WUE; determined from 
images and traditional IRGA measurements (Fig. 9) illus- 
trates the robustness of this imaging approach for rapidly 
assessing WUE,. Although not statistically significant, the 
imaged values in Fig. 9 tended to be in general lower than 
those measured directly. There are two plausible explana- 
tions for this; first, measurements were not taken on identical 
plants and, secondly, and more probably, these small differ- 
ences are the result of a higher boundary layer conductance 
in the IRGA compared with the imaging system. 

Although the potential to couple chlorophyll fluorescence 
and thermal imaging techniques has been explored in several 
studies (Omasa and Takayama, 2003; Chaerle et al, 2007, 
2009), most of these have been conducted at the small scale 
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Fig. 9. A comparison of IRGA measurements of WUE, and 
calculated values of /WUE, from images captured using the 
combined imaging system. WUE, was measured from leaves 
(circles) of A thaliana during a step-wise change in light at a C0 2 
concentration of 400 jimol mol" 1 (see Fig. 8); data are means with 
standard errors (n=5-8). Measurements were taken from individual 
leaves (squares) at C0 2 concentrations between 1 00 umol mol" 1 
and 2000 jimol mol -1 . Air temperature and VPD were 25 °C and 
1 .2 kPa, respectively. The solid line represents a 1 :1 relationship. 

(Omasa and Takayama, 2003; Messinger et al. , 2006) and have 
focused on evaluating stomatal behaviour relative to photo- 
synthetic performance (e.g. Chaerle et al. , 2005, 2007), includ- 
ing stomatal patchiness (Terashima et al, 1988; Terashima, 
1992; Beyschlag and Eckstein, 1998; West et al., 2005), pho- 
toinhibition treatment (Badger et al., 2009), and abiotic and 
biotic stress (Omasa and Takayama, 2003; Aldea et al., 2006). 
Morison et al. (2008) were the first to highlight that the F q '/F m ' 
ratio from fluorescence, along with an indicator proportional 
to g s (7 g ; see Jones, 1992) from thermography, provided the 
possibility to screen plant material non-destructively for A 
and the assimilation transpiration ratio (ATR; a parameter 
mathematically equivalent to WUE). Several previous stud- 
ies have overlaid images of chlorophyll fluorescence and leaf 
temperature to detail the relationships between A and g s . 
However, crucially, this is the first study that has converted 
image data to values of A and g s (using well-defined calibra- 
tions) and used these data to produce quantitative images of 
WUE; collectively. 

An important potential application of the imaging system 
is to screen plants for differences in WUEj. A demonstra- 
tion of this potential is shown in Fig. 6 where clear differ- 
ences in /WUEj images between three WT A. thaliana and 
OST1-1 mutant plants were observed. Images of A and g s 
from these plants indicate that differences in g s , rather than 
A, between the WT and the mutant are primarily responsi- 
ble for the differences in /WUE;. As WUEj is a function of 



both photosynthesis and stomatal behaviour, it is essential 
that future screening and selection of plants with improved 
WUEj is not at the expense of overall carbon gain that 
may translate into reduced crop yield. The ability to deter- 
mine whether differences in A or g s account for differences 
in /WUEj is essential for understanding the physiological 
mechanisms limiting WUEj, as high values of WUEj can also 
be achieved with low A and g s , highlighting the importance 
of developing a system that can measure these two param- 
eters independently. The effects of the differential contribu- 
tion of A and g s to WUE; is exemplified in Fig. 6 showing 
a higher stomatal conductance in the OST mutants com- 
pared with the WT at 20min, but with no net gain in net 
C0 2 assimilation, which resulted in a reduced WUE in these 
plants. Additionally, after 35min, g s had increased in both 
the WT and mutant plants, although this increase was sig- 
nificantly greater in the WT plants, resulting in a greater A. 
The corresponding images of 7WUE; at this time point are 
in general lower than those observed 1 5 min earlier, and the 
initial advantage observed in the WT plants has been lost 
although an overall greater carbon assimilation rate is appar- 
ent. These data provide a prime example of the importance 
of assessing the phenotypic components that drive WUE; 
in screening approaches and protocols to ensure the appro- 
priate combination of physiological traits is selected for 
improved WUE;. The ability to quantify A from images of 
F q '/F m ' is only possible due the built in capability of switch- 
ing from atmospheric to low 0 2 in the measuring chamber 
which allows a linear relationship between F q '/F m ' and A to 
be observed, facilitating rapid comparative measurements 
of A, g s , and /WUEj. Six plants were imaged by the system 
in Fig. 6; however, it would be possible to image a greater 
number of smaller plants but with reduced pixel resolution. 
Lowering the 0 2 concentration to 20 mmol mok 1 has little 
effect on g s over the short time periods required for the meas- 
urement of F q '/F m ', with stomatal responses only apparent 
after ~5-10min. However, it should be noted that determin- 
ing WUE under non-photorespiratory conditions may differ 
if, for example, mesophyll conductance was different in one 
specimen relative to another under investigation. 

Generally protocols for screening differences in A and g s 
are made at steady state (Merlot et al., 2002; Schurr et al, 
2006; Fiorani and Schurr, 2013); however the potential for 
plants to modify A and/or g s in changing environments, 
such as those found in the field, may well be important in 
determining optimal productivity. The imaging system was 
designed to facilitate rapid changes in C0 2 and 0 2 concentra- 
tions, humidity, PPFD, and temperature to generate dynamic 
responses of A, g s and /WUEj (Figs 7, 8). The success of this 
type of dynamic screening protocol depends entirely upon 
the ability to control the measurement conditions tightly, 
which does not usually take priority in many of the existing 
phenotypic platforms. 

In conclusion, this novel imaging system provides the 
opportunity rapidly to assess spatial and dynamic differ- 
ences in A, g s and WUE; in multiple plants under well-defined 
environmental conditions. This should facilitate improve- 
ments in the throughput of plants in phenotyping and 
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screening protocols and consequently in the development of 
programmes for improved crop productivity. 

Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Pixel value distributions for images of C0 2 
assimilation (a and d), stomatal conductance (b and e), 
and 7WUE, (c and f) for a single WT (black) and OST 
mutant (grey) at 20min (T 20 ) and at 35min (T 35 ) under 
200 umol nr 2 s" 1 PPFD (see also Fig. 6). 

Figure S2. Pixel value distributions for images of C0 2 
assimilation (a), stomatal conductance (b), and IW\JE 1 (c) for 
a single plant during a step-wise increase in PPFD (see also 
Figs 8 and 9). 
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